Wave-packet view of the Dirac electron 
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By viewing the electron as a wavepacket in the positive energy spectrum of the Dirac equation, 
we are able to achieve a much clearer understanding of its behavior under weak electromagnetic 
fields. The intrinsic spin magnetic moment is found to be established from the self-rotation of the 
wavepacket. A non-canonical structure is also exhibited in the equations of motion due to non- 
Abelian geometric phases of the Dirac spinors. The wavepacket energy can be expressed simply 
in terms of the kinetic, electrostatic, and Zeeman terms only. This can be transformed into an 
effective quantum Hamiltonian by a novel scheme, and reproduces the Pauli Hamiltonian with all- 
order relativistic corrections. 
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When the electron spin was first discovered, Uhlen- 
beck and Goudsmit regarded that it came from the self- 
rotation of the electron charge sphere 1 -. However, this 
idea was soon swept off by Lorentz, who argued that the 
surface of sphere would rotate with a tangential speed at 
137 times of the speed of light in order to produce the 
spin angular moment of ft/2 2 .. One is thus forced to ac- 
cept the spin as a discrete degree of freedom in abstract 
form, and eventually finds comfort in the Dirac theory^, 
which explains the atomic spectra with great accuracy. 

In the non-relativistic limit, the electron spin manifests 
through its physical attribute, an intrinsic magnetic mo- 
ment, called the Bohr magneton. It is responsible for the 
Zeeman splitting of atomic spectral lines in a magnetic 
field. The spin also couples with orbital motion in the 
presence of an electric field, such as the Coulomb field in- 
side an atom, which is responsible for the fine structures 
of the atomic spectra. Such effects are nicely captured in 
the Pauli Hamiltonian, which can be rigorously derived 
from the Dirac equation through the Foldy-Wouthuysen 
(FW) transformation 4 . However, the appearance of the 
intrinsic magnetic moment for the spin has been mysteri- 
ous, and the usual physical picture linking the spin-orbit 
energy to the Zeeman effect has been rather cumbersome: 
A simple reference- frame argument misses a factor of two, 
and Thomas precession has to be taken into account to 
make the story complete^. 

In this paper, we take an alternative point of view to 
regard the electron as a wavepacket in the positive energy 
spectrum of the Dirac equation. The wavepacket has a 
minimum radius of the Compton wavelength, which is 
larger by a factor of 137 than the so-called classical ra- 
dius of the electron, the size that Lorentz used in his 
argument. In our viewpoint, the wavepacket is found 
to be self-rotating, with an orbital angular momentum 
twice of the spin, and produces a magnetic moment of 
the Bohr magneton with its circulating charge. There- 
fore, the heuristic picture of Uhlenbeck and Goudsmit is 
essentially correct, and the so-called intrinsic magnetic 
moment of the spin originated from the orbital motion 
at the fundamental level. 



In the presence of weak electromagnetic fields, it is pos- 
sible to formulate a semiclassical dynamics for the vari- 
ables of the center of charge, kinetic momentum and the 
average spin&i'^. The magnetic moment gives the Zee- 
man term as a correction to the wavepacket energy, but 
the spin-orbit energy is notably absent. The spin-orbit 
term can arise only because of a shift from the physical 
position variable to an effective but unphysical variable 
in the Pauli Hamiltonian&i^, which is closely related to 
the non-abelian geometric phase of the Dirac spinors^ 1 -. 
In this work, we formulate a generalization of the Peierls 
substitution 1 - 2 , which transforms the wavepacket energy 
to the Pauli Hamiltonian for arbitrary momenta— by a 
novel quantization procedure. In our treatment the spin- 
orbit coupling is a natural product in the canonicaliza- 
tion of the physical variables, and its form is found to 
be gauge dependent. We also suggest that, for a parti- 
cle with anomalous magnetic moment, its mean position 
calculated from the expectation value of the position op- 
erator is actually not the true center of charge, and a 
correction term is required in predicting the correct tra- 
jectory of the particle. 

We now proceed with the calculations leading to the 
above results. We first construct a wavepacket from the 
positive part of the energy spectru m 14 i 15 , 
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where rji are the probability amplitudes of finding the 
electron in the two bands, and the Dirac electron 
wavepacket constitutes of linearly combination of spin 
up and spin down amplitudes of a free particle with a lo- 
calized distribution function a(q, t) centered at q c in mo- 
mentum space. The phase of a(q, t) is chosen to ensure 
that the mean position of the wavepacket is located at r e . 
The conventional Dirac spinors |itj) are obtained from 
a gauge rotation e' 30 ^/ 2 from (1000) T and (0100) T , 
in which a and (3 are the familiar Dirac matrices, and 
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The intrinsic minimum size of electron wavepacket can 
be obtained by introducing the projection operators, P 
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FIG. 1: (Color online) Distributions of (a) particle den- 
sity p = iw'(r)iu(r) and (b) particle current density j(r) = 
w* (r)caw(r) on the X-Y plane for an electron wavepacket, 
the spin is oriented out of the plane. The units are in A c . The 
colors indicate the magnitude of the density from low (blue) 
to high (red). The (c) and (d) are profiles of (a) and (b) along 
the X-axis. 



and Q, for the electron and positron parts of the energy 
spectrum respectively, with P + Q = 1 . The mean square 
radius of the wavepacket can then be written as^ 



(Ar) 2 = (rPr) - (r) 2 + (rQr), 



(2) 



the wavefunction \w) is an eigenstate of P with eigen- 
value 1. Using this fact, one can see that (rPr) — (r) 2 
is simply the mean squared radius for the projected po- 
sition operator PrP, and therefore is positive definite. 
The remaining part 



(rQr) = J d 3 q\a(q) 
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is also positive definite, where the Co mpton wavelen gth 
A c = h/mc is involved and e(q) = + (%q/mc) 2 . If 
we consider a wavepacket with a relatively sharp shape 
in momentum space, then (rQr) 1 / 2 = 2 <${q ) > wn i cn re ~ 
duces to A c /2 as q c — ► 0. This may be regarded as 
the intrinsic size of the non-relativistic electron. If one 
wishes to probe length scales smaller than the Compton 
wavelength, the positron spectrum must be involved in 
the wavepacket and therefore the physical entity being 
probed is no longer a pure electron. 

Within this framework, the electron wavepacket is 
found to be rotating about its center of charge in gen- 
eral, and this rotation is responsible for the intrinsic mag- 
netic moment usually associated with the spin. The Fig.l 
shows the probability distribution and the current den- 
sity of an electron wavepacket with a spin oriented out 



of the plane in the frame of wavepacket. An orbital cir- 
culation around the spin axis can be clearly seen from 
this figure. The mechanical angular momentum for the 
orbital circulation can be calculated from the formula 
L = m((r — r c ) x v) = m J d 3 r(r — r c ) x j(r). A gen- 
eral analytic formula for this quantity has been derived 
in Ref. 14; , which has the following form L = Cr/, where 
7/ = (rji,rj2) T is a two-component spinor, and 

d 

£ y(q c ) = Yl Muili-Q^ln) x (ui\ca\ Uj ). (4) 
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Notice that indices i and j refer to positive energy levels, 
where index I runs over the negative energy levels. A 
straightforward calculation of Eq. (|4]) yields 



C = 



cr + X 



2 q c 
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(5) 



where ar are the Pauli matrices. It can be shown that a 
natural connection exists between the self-rotating angu- 
lar momentum and the projected spin, 



^ T y(q c ) = h(Ui\H\uj) = eCij(q c 



(6) 



where S are 4x4 spin matrices. At extremely high 
velocity(e — » A c q c ), C reduces to zero but r approaches 
the projection along the direction of momentum (q c -cr)q c . 

The magnetic moment of the wavepacket is similarly 
defined as M = — |((r— r c ) x v) = — | J d 3 r(r — r c ) xj(r). 
This differs from the mechanical angular momentum of 
self- rotation by the factor of e/2 in place of m, which 
ensures the usual relation M = — (e/2m)L to be hold. 
One thus has 
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where g = 2. It is very interesting to see that the so- 
called intrinsic magnetic moment of the spin really comes 
from the self-rotation of the wavepacket. In this sense, 
Uhlenbeck and Goudsmit's view of the electron spin as 
a rotating sphere is meaningful. The sphere is rigid near 
the core in the sense that the rotation velocity of the 
wavepacket is approximately proportional to the radius 
near the center. 

The spin magnetic moment also appears naturally in 
the energy of the wavepacket when a magnetic field is 
present. In the Dirac Hamiltonian, the magnetic field 
enters only through the vector potential, and there is 
not a Zeeman energy term. By expanding the vector 
potential about the center of the wavepacket, we find 
the wavepacket energy up to first order in the gradients 
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£(r c ,k c ) = E (k c ) - ecj){v c ) + L(k c ) • B 

Zmc 



(8) 



where k c = q c + |A(r c ), A is the vector potential, and 
Eo(k c ) = e(k c )mc 2 . The Zeeman energy in Eq. ([5]) is just 
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the gradient correction in the wavepacket energy. It is 
strikingly observed that there does not exist a spin-orbit 
coupling in the wavepacket energy, which one would nor- 
mally think to be present to first order in the electric 
field. It turns out that the semiclassical dynamics of 
the Dirac electron is non-canonical, and spin-orbit cou- 
pling is deeply rooted in the Berry-curvature field which 
specifies this non-canonical structure. By requiring the 
wavepacket to satisfy the time-dependent Dirac equation, 
one can derive the effective Lagrangiani^, 

L eff =inr? t ^+nk c -R+nk c -f c -^A-f c -^(r c ,k c ), (9) 

where R is the Berry connection based on the gauge 
specified below Eq. (1), R = 2e (e+i) ^ c x ( <T )> m wn i cn 
(er) = rfcr). The equations of motion for the center of 
charge position and momentum, correct to linear order 
in fields, then follow as 



nk c — — eE x B, 
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(10) 
(11) 



the Berry curvature F 
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These equations need to be solved in conjunction with 
the equation for spin precession, 
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(e + l)mc 



X (cr), 



(12) 



which agrees with the BMT equation^. The equations of 
motion, Eqs. (fIU |) .([TTj ) . and (|T2j). are all invariant under 
SU(2) gauge transformation. 

The same equations of motion have been obtained 
from a formal semiclassical expansion to first order in h e , 
which is not surprising because our dimensionless weak- 
field parameters are eEX c /mc 2 and eBX c /mc 2 , which 
are both proportional to the Planck constant. Our re- 
sults are not limited to low momenta and become even 
more accurate for high momenta where the energy gap 
is greater. For example, one can obtain the semiclassi- 
cal but relativistic cyclotron frequency uj c when the elec- 
tron is confined in a plane perpendicular to the mag- 
netic field (say E = 0), lu c — eBcj E(k c ), where E(k c ) = 
Eo(k c ) + e/(2mc) L B is the total energy of the electron. 

It is evident from the equations of motion that the 
position and momentum of the wavepacket do not form a 
canonical pair, due to the presence of the gauge potentials 
R and A. Quantization of the non-canonical equations 
of motion for the semiclassical Dirac electron presents 
an interesting and important physics problem&2£. If one 
can find new variables r and p, such that the effective 
Lagrangian is of the following form, 



ihr) 



dt 



p • f - E(r,p), 



then r and p would naturally be a classical canonical 
pair. To linear order of the fields (and up to total time 
derivatives), this is indeed possible by the transforma- 
tion, 



r c -R(k c )-G(k c ); 



p = Kk c A(r c ) 



2c 



B x R(k c ), 



(14) 



where G a = l/2(dH/dk a ) • (R x B). Conversely, one can 
write 



r c = r + R(tt) + G(tt) 



TV + -B X R(7f) 



(15) 



where tv = p + -A(r). This is entirely analogous to the 
Peierls substitution for the momentum variable. Naively, 
one expects the familiar presence of the gauge potentials, 
R and A, to the position and momentum. However, this 
is not enough in this semiclassical formulation and extra 
corrections are required. The G-term would further shift 
the position operator, but does not influence the velocity 
(to linear order in fields); similarly, the B x R-term would 
shift the momentum, but does not alter the force. 

One can re-quantize the semiclassical Dirac energy by 
promoting r and p in Eq {8} to be quantum conjugate 
variables, and (cr) in the same equation to be Pauli matri- 
ces. The result turns out to be precisely the (relativistic) 
Pauli Hamiltonian accurate to all orders of the velocity 
(/ib is the Bohr magneton 



H(r,p) — e(n)mc 2 — e<fi(r)- 



e(e + l)mc 



That is, one can obtain an effective quantum Hamil- 
tonian by obtaining the semiclassical energy first, fol- 
lowed by using the (generalized) Peierls substituion for 
re-quantization. We emphasize that one would fail to 
reproduce the correct Pauli Hamiltonian if the extra cor- 
rections in Eq. (|15p were not included. This alternative 
approach is simpler and more intuitive when compared 
to formal procedures of block-diagonalization, such as the 
Foldy-Wouthuysen transformation. 

In Eq. (|16[) . the spin-orbit coupling emerges from the 
first-order gradient expansion of the scalar potential, 
dcf>/dr ■ R. This reveals a deep connection between the 
spin-orbit interaction and the non-canonical structure of 
the semiclassical Dirac theory. It must be remembered 
that the position variable r in the Pauli Hamiltonian docs 
not correspond to the true position r c of the Dirac elec- 
tron. The former is sometimes called the mean position, 
but it is not really the mean position of the wave packet. 
In the literature, — eR has often been called an electric 
dipolei^, which couples to the electric field to give rise 
to the spin-orbit energy^. This is unfortunately artifi- 
cial, because its existence depends on an unphysical as- 
signment of the electron position, which depends on the 
choice of the SU(2) gauge. Indeed, the equations of mo- 
tion based on the Pauli Hamiltonian is consistent with 

(13) 

the Dirac theory if and only if one recognizes this fact. 
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If the Dirac particle has an anomalous magnetic mo- 
ment ( AMM) , which may be originated from the coupling 
with quantized electromagnetic field, or from the internal 
motion of the constituent quarks in a nucleon, then in ad- 
dition to the energy in Eq. ([8]) , the semiclassical particle 
acquires an extra energy (a = g/2 — 1), 
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= a^B 
(17) 
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(k e • (<r))k c 
A c " 2 e(e + 1) 



to be ro + Rd, where rrj is the usual position opera- 
tor, and Rq = (a/xe / 'e)(—i/3a) is the displacement re- 
quired to generate the electric field term in Hamm- The 
displacement (Berry connection) is originated from the 
projection of a larger state space (for example, the one 
that includes the quantized electromagnetic field) to the 
Hilbert space of the relativistic particle. Thus, the true 
velocity operator of the particle should be, 



B 



Unlike the energy in Eq. l[8]). the spin-orbit coupling is 
now explicit from the very beginning. It indicates that, 
when the particle has a AMM, the true center of charge 
r^, has been further displaced from r c . In order to ob- 
tain the spin-orbit interaction in Eq. (fl7|) . the amount 
of displacement has to be R' = a(Ajl/2e) k c x (er). In- 
deed, if one replaces the r c and e/(2mc) in Eq. (JSj> with 
r' c = r c + R' and eg/(Amc), then the AMM energy of 
the wavepacket can be fully reproduced from the sim- 
pler semiclassical energy in Eq. ijHJ). Furthermore, fol- 
lowing the same re-quantization scheme using the gener- 
alized Peierls substitution, one can indeed obtain the (rel- 
ativistic) Pauli Hamiltonian with the exact AMM terms 
(Eq. (33) in Ref3). 

We emphasize that, since one expects no electric dipole 
from a point particle, the presence of the electric dipole 
energy indicates that the theory is only an effective one, 
and the position of the particle in the effective theory 
should differ from the true position by a Berry con- 
nection. Of course, in predicting the trajectory of the 
point particle, one always has to refer to its true po- 
sition. Therefore, in the Dirac theory with J?amm = 
afiB(3(— ioL ■ E + X ■ B), the true position r is conjectured 



r D + Rz 
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a/3- 



a 2 ^(ExS 

mc 



7sB). 
(18) 

The magnitude of the anomaly a for an electron is of the 
order of 10 -3 and might be too small for such a difference 
to be observed. However, the anomaly for a nucleon is of 
order one. Therefore, it is possible to observe and verify 
such a deviation in experiments using relativistic proton 
beams 

In summary, we demonstrate an alternative method 
to investigate the Dirac electron by semiclassical 
wavepacket approach which provides a much more in- 
tuitive and pedagogic picture of the Dirac electron. In 
this framework, the intrinsic size of the non-relativistic 
electron is of the order of A c , and the self-rotation motion 
generates a magnetic moment equals precisely the con- 
ventional spin magnetic moment. The dynamics of the 
wavepacket is also studied and reveals a direct link be- 
tween the Berry connection and the spin-orbit coupling. 
By re-quantizing the semiclassical Dirac theory through 
a generalized Peierls substitution, we recover the rela- 
tivistic Pauli Hamiltonian to all orders of velocity. 
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However, when the velocity is projected to the subspace for 
an electron with positive energy, one has PcolP = n/(em). 
It would cancel with the term P[3ir/mP — ir/m to lead- 



ing order of the velocity. Therefore, the leading (field-free) 
correction is of the order of a(v/c) 2 . 



